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Abstract 

Membrane tubes (spicules) arise in cells, or artificial membranes, in the nonlinear deformation 
regime due to, e.g. the growth of microtubules, actin filaments or sickle hemoglobin fibers towards 
a membrane. We calculate the axial force exerted by the cylindrical membrane tube, and its 
average radius, by taking into account steric interactions between the fluctuating membrane and 
the enclosed rod. The force required to confine a fluctuating membrane near the surface of the 
enclosed rod diverges as the separation approaches zero. This results in a smooth crossover of 
the axial force between a square root and a linear dependence on the membrane tension as the 
tension increases and the tube radius shrinks. This crossover can occur at the most physiologically 
relevant membrane tensions. Our work may be important in (i) interpreting experiments in which 
axial force is related to the tube radius or membrane tension (ii) dynamical theories for biopolymer 
growth in narrow tubes where these fluctuation effects control the tube radius. 
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There has been much recent interest in the formation of tubular membrane tethers or 
spicules These arise from the action of localised forces that 

act normal to the membrane and give rise to narrow membrane tubes in the nonlinear 
deformation regime. Such forces can arise from the polymerization of fibers, including 
actin^(| , tubulin 11], or sickle hemoglobin jl^ . into stiff fibers or bundles of fibers. Similar 
structures appear on cell membranes as filopodia or smaller tubular excursions lri | 
including the spicules of stellated or sickled red blood cells containing sickle hemoglobin 
(HbS) 12] or neural growth cones as well as on vesicles observed in vitro 9]. When the 
tube is long (many times its radius) it is approximately cylindrical on average except very 
close to its endsjlOj], see Fig 1. Throughout we will assume that any growth of the fiber and 
tube is quasi-statically slow. 

In the following section we outline a self consistent mean field analysis of the radial 
membrane fluctuations. In this the average radial extent of the fluctuations are controlled 
by the presence of the enclosed rod, an approach that is analogous to Helfrich's highly 

,15]. We will work in units in which k B T = 1. 

In order to describe our cylindrical membrane, we use the Hamiltonian H = He + H$ 
with 

H E = f[a + f c 2 ] y/g d(f)dz - fL ^ 



H s = J[A + J< 



z ) ~ r ) + 5 ( r ( 



-fL 

z) — r) 2 ] d&dz 



l4 lid Il7 |. containing both 



where He is the usual Hamiltonian for membrane elasticity 
surface tension (a) and rigidity (k) controlled terms [18], the latter varying with the square 
of the local membrane curvature c. We have also included in He an axial force term, /, 
which arises from the polymerizing fiber, and controls the axial length, L, of our membrane 
tube. In what follows we neglect any effects that could arise from leaflet asymmetry or 
spontaneous curvature. H$ contains a harmonic potential (with strength ~ C) that confines 
the size of the membrane fluctuations about the average tube radius f. It also contains a 
'radial force' (or pressure) term (with strength ~ J), which controls the average radius f 
of the membrane tube and arises from the asymmetric nature of the constraint provided 
by the rod. This term can also be used to include any hydrostatic or osmotic pressure 
differences between the inside and outside of the tube although we neglect these in what 
follows. Additionally, H$ contains a term involving A which is convenient for normalisation 
of the steric potential. This, most general harmonic potential, will be used to model the steric 



2 



interactions between our membrane and polymer rod by way of a mean-field approach. An 
analogous treatment has proven to be remarkably successful in describing the steric repulsion 



between flat membranes 



14. 



15(. We proceed from Eq (J) by writing r(0, z) = r + £r(0, z) 
and expanding the energy H to quadratic order 

H Q, Q 

in the radial perturbation 

5r((j), z) about the average tube radius f. A convenient Fourier representation of the radial 
fluctuations is 

6r(<f>,z)= ^2 ^2 fr ™ ex Pl™^ — ] ( 2 ) 

n=— oo m=~ oo 



With q = ^r 2 - ) this yields the energy which we write as a perturbative expansion H 
H + 5H + 5 2 H + ... with 

H Q = — + 2iraLf + 2ttAL - fL 

r 

SH = =^< 5r>=0^>J = -a + 

2r l 

1 2nnL " 



$ H — _ 3 ^ ^ ^ ^ ^ r nm^ r —n—mK nm (3) 



n=— oo m=— oo 



involving a kernel 



i^ m = (™V + ™ 2 ) 2 -^(l-^) 

5m 2 / i f 2 | i 1 _i_ Cf 3 



(4) 



which, in turn, involves the radial length scale r = v/^- The first order perturbative 
contribution is required to vanish so that f indeed represents the true average (or ground 
state) membrane tube radius j^]. This condition then implies (from Eq (J3J)) that the 'radial 
force' is J = —o + Sx. The quadratic fluctuations in the radial displacement (around f), 
contribute at order 5 2 F, and depend on the strength of the harmonic potential in Eq ((TJ, 
via the parameter C. The presence of the rod sterically constrains the membrane radius, 
r((p,z), to remain always greater than the rod radius b, see Fig 1. It thus follows that the 
presence of the rod has an effect on both the average radius, f, and the fluctuations, 5r((j), z). 
The mean squared amplitude of the fluctuations, (Sr 2 ), depends on the parameter C, as can 
be seen from Eq (J3J). We must now determine this self consistently. Finally it may help to 
note that the rod, and the potential that we employ to mimic it, have the additional role of 
suppressing instabilities t 
tubes (such as 'pearling' 



rat are known to occur on certain 'rodless' cylindrical membrane 



24 |2a |2£( for example) 



The free energy of the tube F = Ho + AF can be shown to involve the perturbative 
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correction 

AF = 

2 



oo oo 

AF = o E E ln ( 5 ) 



n=— oo m=— oo 
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one must 



However, as is typically the case in functional integrals of the above kind 
invoke a process of 'renormalisation' in order to remove all divergent (small-wavelength) 
contributions to AF, and hence the free energy F. This process is typically achieved via 
appropriate 'counter-terms' in F that subtract from AF those parts (and only those parts) 
which are divergent, while leaving behind a physically meaningful finite contribution 
Moreover, just such a candidate 'counter-term' exists in our calculation of F, namely the 
parameter A, which arose in the most general steric potential Hs- Furthermore, one can 
physically motivate an explicit choice for A as follows. We aim to calculate and compare the 
free energy difference between the case when the rod is present and when the rod is absent. 
When the rod is absent the steric harmonic potential (of strength C) vanishes, as do terms 
involving C that appear in K nm . So we choose the parameter A so that in the limit C — > 
we obtain AF — > for consistency. We must then choose 



A 



^ oo oo 

— E ^(K nm \ c=0 ) (6) 



A-kL 

n=— oo m=— oo 



After we have integrated out all radial membrane fluctuations we therefore obtain 

F ^ + 2mL ^ fL + l± f>(^) (T, 

n=-oo m=-oo x ' ' 

Further, by converting the above summations into integrals, and performing a change of 
integration variables, we can write our final expression for the free energy as 

F = + 2naLf - fL (8) 

r 

in which we have defined an 'effective' membrane bending modulus through which one may 
interpret all of the effects of the steric potential 

J /"27T /*00 

K eff = K + 7T^ d6 dp\og(l + Cf 3 /(KP)) (9) 
71 " Jo Jo 

with 

P = p 2 -^l- T ^J -2pcos 2 6 + l (10) 

We can now state quantitatively the physical condition that we wish to impose on our 
membrane to mimic the steric influence of the rod 



This gives the necessary self-consistency condition for the strength of the harmonic potential 
given that 

_3 oo oo 

= ^-T E E *™ ( 12 ) 



n=— oo m=— oo 

By again converting the summations into integrals and changing variables, we arrive at the 
following consistency equation for C 

b\ 2 i r 27r r°° 1 

In order to calculate the average radius f for the membrane tube we merely need to minimise 
F by setting ^ = 0. Similarly the force / required to maintain the axial length L of our 
membrane tube follows from ff- = 0. 



Wide tubes ( f/b » 1 , Cr 3 /K < 1 ) 

This regime corresponds to the case when the equilibrium radius of the membrane tube 
is typically much larger than the radius of the enclosed rod (f/b 3> 1)- The steric effects of 
the rod in this limit should therefore be weak (Cf 3 /n <C 1), and as a consistency check, we 
recover the results of jlj], namely f = r = and / = Ih^Iok. For the present purposes 



2<T 

^2 



we can expand about r = tq using Eq ()13j) to obtain C = ^pj- exp (— 47r/tv / 2) • A similar 
exponential form for the confining potential is found for stacks of flat membranes at large 



inter-membrane separation 
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3l| . Substituting into Eq we obtain to leading order 



2vr 2 / \ 
K e ff/K — 1 + exp ^-47tkv / 2J (14) 

From which we can see that the rod/membrane steric interaction produces a small correction 
to k in this weak confinement limit (b/r <C 1). The energy F is given by Eq (jHJ) giving the 
axial force / = F/L and, by minimisation, f to leading order 

f/r = 1 + y exp (-47^^) (15) 

recovering the well known limiting results [1], plus the leading order steric correction term. 

The axial force is then / = 2-k\/2ok (l + ^-exp (— Attk\/2) j , which corresponds to the 
well known result Q but includes the leading order correction term due to steric effects. 
Given a typical value of k ~ 10 (in units of k B T), we can see that the steric correction terms 
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in f and / become very small in the present limit (6/ r o 1)- However, as b/ro — > 1, we 
cross over to the narrow tube regime in which there is strong confinement of the membrane, 
as discussed below. 



Narrow tubes ( r/b ~ 1 , Cr 3 /n 3> 1 ) 

This case corresponds to when the average radius of the membrane tube is almost equal 
to the radius of the enclosed rod {r/b ~ 1). The steric effects of the rod in this limit 
should therefore be very strong, and the strength of the self-consistent, confining, harmonic 
potential becomes very large (Cf 3 /^ ^> 1). In this case we can approximate Eq (fT3*j) as 
follows : 



(f-b) 2 i r 2n 
1 ' - 1 de I dp 



b 2 SttWo Jo V + ? 

as C — > oo (16) 



8 V nCb 3 

Hence in this limit C = ^ ^-b) A • Substituting this value of C into Eq @ for « e //, we obtain 
to leading order (as C — > oo and f — > b) 

K ^ ,K = l + ^ 2 j¥hf (17) 

From Eq ()17jl we can see that the rod/membrane steric interaction can dominate the effective 
rigidity in this strong confinement limit. Furthermore, substituting this dominant value of 
K e ff into Eq (JSJ) we obtain (to leading order) 

F = 2naLb + 2naL(f -b)-fL+ ^J~^~ 2 (18) 

A contribution to the energy that scales as the inverse squared distance, similar to the 
one appearing here, is well known for flat, parallel membranes at small inter-membrane 



separation 
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Proceeding as before we find r = b + (3^) 3 , and 



, 3tt (o 2 b\ 

f = 2*ab+- — (19) 
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This, and particularly the ~ a 2 / 3 correction, is one primary experimentally testable predic- 
tion of our study. These results are therefore consistent at high enough tension or rod radius 



nab 2 > 1. The result / « 2nab arises from the work done against surface tension in laying 
down additional membrane area near the surface of a rod of fixed radius b. 

DISCUSSION 

Typical experimental parameter values for biological membranes might be K ~ 

10~ 19 J, a « 10~ 4 Jm~ 2 which give rise to a typical tube radius of r ~ 20nm. Comparing 
this value for the membrane tube radius with typical biological fibers (or fiber bundles) with 
radii b ~ 10 — lOOnm, we conclude that steric effects will be important and our results 
are therefore likely to be of biological significance under physiologically relevant conditions. 
Typical forces are also found to be in the experimentally relevant range of / ~ 10 — lOOpN. 
The steric corrections to f and / when the membrane is strongly confined by the presence of 
the enclosed rod are governed by the dimensionless parameter ~ {a Kb 2 ) 1 which is significant 
for the typical values given above. Hence additional steric corrections may be expected to 
contribute a measurable correction to the expected radius f and force /. Our work may be 
of significant importance since measurements of the axial force are often used as an indirect 
way of estimating the tube radius. 

CONCLUSIONS 

In this work we have studied the effect of an enclosed rod on the fluctuations of a cylin- 
drical membrane tube. Such tubes, or spicules, have been observed in various biological and 
chemical experimental systems. We have calculated the axial force exerted by the tube and 
its radius via a self-consistent harmonic potential that models the steric constraint of the 
rod on the membrane fluctuations. The axial force diverges as the tube radius decreases and 
the membrane approaches the surface of the enclosed rod. Our approach is similar to to the 
approach used for flat membranes Q|. However, in the present case, the radius of the 
fluctuating membrane tube is only constrained from inside by the presence of the rod and 
not from outside. This lack of symmetry gives rise to two parameters (C and J) instead of 
one. We believe that our results will be of importance in the context of polymerising fibers 



within cell membranes 



aaiaa ;34j. In particular we give the first quantitative treatment 



of the crossover between the 'free' membrane tube, with radius r ^> b, for which the axial 
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force is 'small' / ~ \JTio and the f — > b regime, in which the radius of the tube is almost 
constant and / ~ a. Our analysis is also likely to be of significant utility in approaching such 
unsolved problems as the transport mechanism of material within the tube, e.g. monomers 
destined to polymerize at the tip of a growing rodj^j]. Such polymerisation processes may be 
expected to depend intimately on the tube radius in the high tension (narrow tube) regime. 
In this sense future dynamical theories for biopolymer growth in tubes (spicules) are likely 
to rely on an accurate treatment of the fluctuation effects that we have analyzed here. 
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Figure caption 

Fig 1 Sketch of an axial slice through a membrane tube (left) and a section of tube (right) 
of average radius f enclosing a cylindrical rod (of fixed radius b). The radial fluctuations are of 



amplitude 6r(z, <fi) and have mean squared extent (Sr 2 ). 



